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Abstract Microtubules, major elements of the cell skeleton are, most of the 
time, well organized in vivo, but they can also show self-organizing behav- 
iors in time and/or space in purified solutions in vitro. Theoretical studies 
and models based on the concepts of collective dynamics in complex systems, 
reaction-diffusion processes and emergent phenomena were proposed to ex- 
plain some of these behaviors. In the particular case of microtubule spatial 
self-organization, it has been advanced that microtubules could behave like 
ants, self-organizing by 'talking to each other' by way of hypothetic (because 
never observed) concentrated chemical trails of tubulin that are expected to 
be released by their disassembling ends. Deterministic models based on this 
idea yielded indeed like- looking spatio-temporal self-organizing behaviors. 

Nevertheless the question remains of whether microscopic tubulin trails 
produced by individual or bundles of several microtubules are intense enough 
to allow microtubule self-organization at a macroscopic level. In the present 
work, by simulating the diffusion of tubulin in microtubule solutions at the 
microscopic scale, we measure the shape and intensity of tubulin trails and 
discuss about the assumption of microtubule self-organization due to the pro- 
duction of chemical trails by disassembling microtubules. We show that the 
tubulin trails produced by individual microtubules or small microtubule arrays 
are very weak and not elongated even at very high reactive rates. 

Although the variations of concentration due to such trails are not sig- 
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nificant compared to natural fluctuations of the concentration of tubuline in 
the chemical environment, the study shows that heterogeneities of biochemi- 
cal composition can form due to microtubule disassembly. They could become 
significant when produced by numerous microtubule ends located in the same 
place. Their possible formation could play a role in certain conditions of 
reaction. In particular, it gives a mesoscopic basis to explain the collective dy- 
namics observed in excitable microtubule solutions showing the propagation 
of concentration waves of microtubules at the millimeter scale, although we 
doubt that individual microtubules or bundles can behave like molecular ants. 

Keywords Microtubules • Self-organization • Tubulin trails • Ant-based 
model • Microscopic scale model • Reaction-diffusion 



1 Introduction 



In this article we aim to characterize by a model of diffusion working at the 
microscopic level the shape and the nature of hypothetic heterogeneities in the 
free tubulin concentration profile that are expected to be produced by disas- 



sembling microtubules in Robert et al ( 1990 ); Tabony et al ( 1999 ); Glade et al 



( 2002 ) ; Tabony ( 2006 1 . Their existence in the form of concentrated regions 



of tubulin-GTP or tubulin-GDP conditi on the validity of certain models of 
microtubule self-organi zatiorp] in vitro (Pirollet et al||l987 Carlier et al||l987 



et al 2002 



et al||1999) based on reaction-diffusion processes (Robert et al 



Mandelkow et al| [T989; Ta bony and Job||1992b[ |Papaseit et al||1999[ |Tabony 



1990 



Glade 



Sept 1 1999 1 . The numerical experiments realized in this work added 



to experimental observations from the literature bring new elements to under- 
stand microtubule self-organization at the microscopic level. 

Here, we want to answer two fundamental questions: (i) the question of the 
formation of chemical trails by disassembling microtubules and (ii) the ques- 
tion of their relative influence during self-organizing processes. Our results are 
not in favor of the hypothesis of an ant-like behavior either for individual mi- 
crotubules or for microtubule arrays at a microscopic level, although we think 
the question is still open. We however confirm the idea that large scale areas of 
different composition in tubulin (tubulin-GTP vs tubulin-GDP) can form and 
explain large scale self-organizing phenomena in microtubule solutions, such 
as the temporal or spatio-temporal oscillations of microtubule concentration 



described in Pirollet et al (1987); Carlier et al (19871; Mandelkow et al (1989). 



We also make a review of several experimental conditions and results that rein- 
force the thesis of a biomechanical scenario for spatial self-organizing solutions 



as proposed in Baulin (20031; Portet et al (20031; Ziebert and Zimmcrmann 



1 In this article, the term self-organization has a different meaning than self-assembly. 
Self-assembly addresses to the process by which an individual microtubule forms sponta- 
neously by assembly of tubulin-GTP subunits ||Wciscnbcrg ]1972| l. Self-assembly is of course 
a kind of self-organization, but here self-organization designs specifically a more macroscopic 
organizational level. It refers to the order that can appear from particular distributions and 
orientations of microtubule populations in solution. 
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(2004 


) ; Liu et al 


(200G): 


Baulin et al 


( 2007|) , instead of a 'social insects'-based 


process as proposed in 


Robert et al 


(1990); 


Tabony et al 


(1999 


); 


Glade et al 



(2002); Tabony (2006). 



1.1 Microtubules 



Microtubules are biological supra-molecular assemblies with a micrometer scale 
(about 30 nm diameter and several micrometers long). They are energy- 
dependent reacting fibers present in most of the living cells, often in large 
amounts (e.g. they represent 10% of the total amount of proteins in neurons). 
In addition to actin micro-filaments and intermediate filaments, they form the 
cell skeleton (cytoskeleton) that gives shapes and biomechanical properties to 
the cells. Actin filaments and microtubules don't only constitute simple static 
structures ; these fibers react, grow or shrink, forming that way dynamical 
structures that self-adapt to the changes of cell states - particularly changes of 
energy levels - and to the mechanical, electrical or magnetic exogenous stimuli 
(|Vassilev et al|1982||Tabony|1994[|Tabony and Job|1992a[|Papaseit et al|2000 



Glade and Tabony||2002f |2005| |Glade et al||2006| ls.J.Crawford-Young 



2006 



Roesner et alp006 GalimbertTeFap006 ; Kroupova et~alp007; Coleman et a 



20071 |Tabony et al||2007| |Ingber|2008| |Wang et al| 2008; Qi an et al|2008| [Yang 
et al 2008| |Sieberer et al||2009| |Qian et al||2009a|b| . Microtubules get also 



involved in numerous cell functions, often constituting organelles : they form 
the centriole, the organizing center of the cell from which the microtubules 
radiate, the mitotic spindle that drives the chromosomal segregation during 
cell division, the elastic 'motorized arms' of cilia and flagella. They also serve 
as rails guiding the molecular motors that drive the active transport of phos- 
pholipidic vesicles. 

Microtubule - as actin filaments do - are able to self-assemble from their 
constituting bricks called tubulin heterodimers (composed of a and /? tubulin) 
( Weisenberg|1972" ) once warmed from 4 ° C to about 35 ° C (this range of tem- 
peratures constitutes an optimum, but microtubules can also form at room 
temperature at slower rates). Tubulin heterodimers are associated either with 
GTP (guanosine triphosphate) or GDP (guanosine diphosphate). Both can as- 
semble but the assembly of tubulin-GTP is more favorable than tubulin-GDP. 
Microtubules are usually extremely reactive and weakly stable : their ends 
assemble or disassemble quasi permanently without periods of stability. Their 
reactivity is, in first approximation, linked to the amount and the nature of the 
tubulin heterodimers encountered locally in the neighborhood of their ends: 
the reaction of assembly is more favorable when the medium is concentrated in 
tubulin-GTP heterodimers than in tubulin-GDP heterodimers. The reactivity 
of microtubules also results from the very complicated intrinsic mechanochem- 
ical dynamics that occur within their ends (VanBuren et al 2005). Due to 



those non-linear dynamics, they can show either stochastic-looking individual 
behaviors - also called dynamic instability ( Mitchison and Kirschner]|1984 ) -, 
where microtubules switch apparently in a stochastic manner between rapid 
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phases of assembly and of disassembly, or pretty regular ones such as the well- 
described treadmilling process (microtubules assemble regularly at one end 
while disassembling at the other, causing that way the microtubules move in 
their direction of growth) ( Margolis and Wilson|1978 ). Moreover, assembly or 
disassembly reactions modify locally (at the microtubule scale) or more glob- 
ally (depending on the ratio between reaction rates and molecular diffusion) 
the respective amounts of the different chemical species in the medium {e.g. 
tubulin-GTP, tubulin-GDP, GTP, GDP, ...). These c hanges can have in re- 
turn an effect on the microtubule population dynamics ( Mandelkow et al|1989 



Jobs et al 1997 Sept 1999 Deymier et al 2005). A synchronizing of micro- 



tubule reactivities in time and/or space can arise from this coupling between 
the two chemical states of tubulin, i.e. assembled tubulin or free tubulin het- 
erodimers, notably under the control of the regeneration rate of tubulin-GDP 
into tubulin-GTP (|Melki et al|[l988]). 



1.2 Behaviors of microtubules in vitro. 



In vitro solutions that only contain purified tubulin and GTP (the energy 
source of the reaction) in a buffer composed of various ionic species (Mg 2+ 
notably) can show singular macroscopic behaviors depending on the reactive 
conditions (temperature, ionic concentrations, tubulin concentration ...) (Ta- 
bic [lj. In such solutions, one can observe behaviors from 'simple' temporal 
periodic oscillations of microtubule concentration within the whole solution 



(Pirollet et al 1987 Carlier et al 1987 Melki et al 1988 Mandelkow et al 



1988 Hitt et al 1990), to spatial 'stationary' morphologies with very com- 
plicated structures (Tabony et al 1999). Intermediate behaviors also exist, 



showing concentration waves of microtubules forming and propagating peri- 
odically throughout the solution ( Mandelkow et ai||1989 ). This self-organized 
spatial behavior looks similar to those observed in some excitable media such 
as the Belousov-Zhabotinskii reaction ( Vavilin et al|1967b|a Zhabotinsky and 
Zaikin||1973| |Zhabotinsk"ii||1974 1 . 

Other varieties of spatial self-organization were reported in similar solu- 
tions in vitro of purified microtubules containing also microtubule associated 
proteins ( Nedelec et al|1997 Surrey et al|200l] ). In the cells, microtubules are 
frequently associated with such proteins (MAPs). Most of them are molecular 
motors {e.g. kinesin, dynein) that use the energy of phosphated nucleotides 
(GTP) to propel and move directionally at the surface of microtubules while 
carrying vesicles. They can also link individual microtubules between each 
others in such a way they form a mechanically coupled network. Some have 
a stabilizing effect on microtubules. To explain how these particular solutions 
self-organize, the authors use numerical models based on the coupling between 
microtubule reactivity and mechanics biased by the reactivity and mechani- 
cal constrains induced by molecular motors ( Surrey et al||200T Nedelec||2002 



Nedelec et al|20 03 ) . Their numerical simulations yield very realistic behaviors 
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suggesting that their model correctly explains how microtubules and motors 
self-organize. Similar works succeed in explaining the spatial organization of 
actin networks (Rolland et al 2008 Pollard and Berro 2009 Reymann et al 



2010) 



The behaviors described by 


Pirollet et al 


(19871 


Carlier et al 


(1987 


). Melki 


et al 


(1988) 


Mandelkow et al ( 


1988 


1989 


)• 


Hitt et al 


(|1990) and|Tabony and 


Job 


1992b) 


; Papaseit et al 


(1999); 


Tabony et al 


(] 


_99< 


i ) however imply differ- 



ent processes because they don't contain any molecular motor (except in some 



experiments in Hitt et al ( 1990)). The cases of temporal oscillations or spatial 
waves of microtubule concentration are not yet ambiguous. They clearly imply 
the major contribution of microtubule dynamics and their coupling through- 
out the solution by diffusion of matter (at least tubulin heterodimers) . Indeed, 
these alternating succession of assembling and disassembling periods were di- 
rectly observed by dark field microscopy at the mesoscopic level (observation 
field of about dozens of microns large, largely greater than microtubule di- 
ameters) ( Mandelkow et ai]|1989 ). Further, several models based on reactions 



coupled to diffusion processes were proposed (Carlier et al 



2002 



1987 



|Robert et al]|i~99"0l |Marx and Mandelkow||1994[ |Jobs et al 



Glade et al 



1997 Sept 



1999 ) and some were able to reproduce the observed behaviors. In the article 
we will refer frequently to 3 different typical situations - (i) temporal oscil- 
lations, (ii) spatio-temporal waves, and (iii) stationary morphologies - first 



described respectively by (i) Pirollet et al| ( 1987 ) and Carlier et al 
Melki et al|(|1988| ) and |Mandelkow et al| ( |1988^ and (iii) |Hitt et at 
Tabony and Jobl (|1992b|). 



19871, (ii) 



( 1990 ) and 



Stationary self-organized microtubule morphologies are the more compli- 
cated to explain. Indeed, they don't clearly show as important dynamic be- 
haviors as those described before, as to say dramatic variations of the con- 
centration of microtubules in time (microtubule concentration is followed by 
turbidity measurements at 350 nm) or rapid wave propagations. In these so- 
lutions, the spatial self-organization is slow and progressive. Morphologies 
develop during several hours (an average of 5 hours is required for obtaining 
well organized morphologies, although the self-organizing process continues). 
After a unique 'overshoot' at the initial stages of the reaction (after about 5 



minutes), the microtubule concentration stabilizes (Tabony 1994 ). The process 
is obviously energy dependent (since tubulin-GTP is necessarily hydrolyzed for 
the microtubules to grow) (Tabony and Job 1990 Tabony 1994) and seems 
to be due to a combination of coupled diffusion and reactions (Papaseit et al 
1999 Tabony et al 1999[ ). The morphologies obtained are very reproducible, 
and complicated : they present alternating stripes of concentration (Papa- 



seit et al 1999) and orientation (Tabony and Job 1992b Tabony 19941 of 



microtubules at a macroscopic level (millimeter) and moreover show nested 
sub-levels of self-similar organization between the macroscopic level and the 
level of microtubule bundles ( Tabony 1994[ ) . Finally, it has been also shown 
that the presence of a weak external field such as gravity (Tabony and Job 



1992a| |Papaseit et al||2000| |Glade and Tabony||2002[), magnetic fields ( |Glade 
and Tabony 2005 Liu et al 20061, or vibrations (Glade et al 2006) at the 
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early stages of the reaction is necessary to trigger the development of these 
self-organized morphologies at a macroscopic level. Living cells and organisms 
are also sensible to weak external fields and this often implies microtubules 
(S.J. Crawford- Young 2006). Coherently, the understanding of the processes 
involved in vitro was though to give information on the manner some external 
fields could act on living cells. Once again, as for microtubule traveling 
waves ( Mandelkow et al|1989 Sept 1999), microtubule stationary spatial self- 
organization has been thought to be an 'emergent' behavior; a consequence 
of the collective dynamics of microtubules over space and time. Moreover, 
the action of the weak external factors on these solutions is understood in 
the sense of a symmetry breaking in a 'complex system' ( Tabony||2006 ) . This 
approach was quite original for the reason that it was an alternative to the 
other plausible advanced explanations of this phenomenon of gravisensitivity. 
Usually, people indeed consider the biomechanics (bundling and buckling of 
the fibers) and/or the various static interactions between microtubule rods 
coupled to growth as a possible way by which they could 'feel' the gravity or 
other external factors and self-organize at a macroscopic level JHitt et al 1990 



Portet et al|2003[[ Ziebcrt and Z immermann|2 004 Baulin 2003; Liu et al 
Ingber||2008| |Baulin et al||2007[ ). 



2006 



1.3 The question of microtubule spatial self-organization. 



The following assumptions are intuitive : the chemical activity of individ- 
ual microtubules or microtubule arrays should cause the formation of local 
variations of concentration and composition of the chemical medium around 
their reacting tips. Such local heterogeneities should affect the reactivity of 
microtubules present in the neighborhood. This could be a way by which mi- 
crotubules 'communicate'. 

It has been proposed repeatedly that the formation of such variations 



could influence microtubule dynamics and self-organization (Tabony et al 1999 
Robert et al||1990| [Glade et al||2002[ ). In 1990, Robert et al published a very 



simple chemotactic model of microtubule self-organization where individual 
microtubules coordinate among themselves and self-organize, following the 



gradients of tubulin concentration self-produced by their own activity ( Robert 



et al 1990). Tabony published another article that reinforces the idea that 
microtubules behave like ants (Fig. [lj, self-organizing by collective dynam- 
ics, with emergent ant-like behaviors such as stigmergy (Tabony 2006). The 



author proposes that this type of physicochemical process, described at the 
microscopic level, is implied in living systems, explaining particularly their 
sensibility to weak external fields or the organizing processes that occur in the 
early stages of embryogenesis. Other recent studies on the effect of magnetic 
and electromagnetic fields or of the weightlessness on cells or microtubule so- 
lutions, works that propose biotechnological solutions based on microtubules, 
but also more general articles on self-organization, refer to these works on 



microtubule self-organization in vitro (S.J. Crawford- Young 2006 Galimberti 
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et al|2006||Yang et al|2006| |Roesner et al|2006] |Kroupova et al|2007[ [Coleman 
et all|2007| |Qian et a.l||2008[ |Coleman et al||2008[ |Yang et al||2008l |Strasak 



et al 2009; Sco ct al 2009; Qian et al 2009a b; Sieberer et al 2009; Qian et al 



20TU1 JTIaiTet al||2010| |01son et al|[20T0l [Johnson and La"ml[20T0l |Moes et al 
2011[ R et al||2011 Ayodele et al||2011[ ). Some still take the reaction-diffusion 



scenario for granted and the others accept the idea that weak external fields 
act similarly on microtubules in vitro and in vivo. 

Unfortunately, the influence of the activities of microtubule ends on other 
microtubules by the intermediate of local variations of the tubulin concen- 
trations - called tubulin trails - remains hypothetic since it has never been 
directly observed experimentally. Moreover, although the concept of reaction- 
diffusion is interesting and certainly significant in some microtubule solutions 



et al 1988 



a macroscopic level (Pirollet et al 


1987 


Carlier et al 


1987 


Melki 


Mandelkow et al 


1988 


1989 


Sept 


1999 


Deymier et al 


2005 


), it 



is probably not reliable in microtubule stationary self-organizing solutions. 
Indeed, in the later, the influence of other effects (e.g. mechanical or elec- 
trostatic interactions between microtubules, molecular agitation at the level 
of individual microtubules) is probably prevalent compared to the effect of 
eventual chemical trails. We do not either think that the macroscopic effects 
observed in microtubule stationary self-organizing solutions under to the ac- 
tion of external fields, can be so easily compared to those occurring in the 
living matter. 

Microtubule self-organization at the level of individual microtubules might 
be based on the existence of local trails or depletions of tubulin concentrations 
generated by the activity of microtubules. These local trails or depletions are 
indeed proposed to be the way by which microtubules 'talk to each other' 
( |Tabony et al|1999[|Glade et al|2002)|Tabony|2006[ ). Here we ask the question 
of how much assembly or disassembly reactions modify the chemical medium 
surrounding the ends of a single microtubule and how does this affects the 
reactivity of neighboring microtubules. 



2 Methods 

2.1 Continuous description of matter and deterministic models of 
microtubule chemical trails. 

A quantitative numerical model of formation of a depleted area of free tubulin- 



GTP around the growing tip of a microtubule was proposed by Odde (19971. 
The aim of this study was to estimate if the formation of this region and its 
homogenization by diffusion could be a limitation for the microtubule growth 
reaction. The result - an analytical solution of a reaction-diffusion equation 
- estimated that for a microtubule growing at 7 /im.min -1 , the concentration 
at the tip is 89% of the concentration far from the tip and that the concen- 
tration gradient is extending to less than 50 nm from the tip (less than 2 
microtubule diameters). The finite differences implementation of partial dif- 
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ferential equation model proposed by Glade et al (2002) also predicted the 
formation of similar depleted areas at the growing ends of microtubules with 
the formation of tubulin concentrated trails at the shrinking ends (Fig. |2j. 
This local phenomenon was hypothesized to be the most fundamental man- 
ner by which microtubules communicate. It was reinforced by observations in 
vitro (Mandel kow et al||1989 ) and in vivo (Keating and Borisy 1999) showing 
that new microtubule bundles preferentially grow in the direction and in place 
of shrinking ones, maintaining that way the existing microtubule bundles. 

Nevertheless, these three numerical models using or describing tubulin dif- 
fusion at microtubule ends ( |Robert et aT||1990| |Odde||1997| |Glade et al||2002[ ) 
were simulating continuous amounts of tubulin, as to say concentrations ex- 
pressed in fimol.l -1 , represented by floating numbers. The use of a continuous 
description of matter (here concentrations) is well-adapted to simulate large 
numbers of reactants. However, it's unadapted to simulate a small number 
of molecules where only discrete descriptions of matter are adapted. In real 
microtubular solutions, only very few tubulin-GDP heterodimers are produced 
by disassembling microtubules and diffuse from their reacting ends, and only 
few tubulin-GTP heterodimers are taken away from the medium and assem- 
bled at the microtubular tips. The use of concentrations in this case is not 
well adapted because, due to continuously-expressed diffusion and reactions, 
it can generate in the simulation tubulin amounts being fractions of individual 
molecules. This is exactly what is observed on Fig. [2j The concentration map 
of tubulin-GTP shows small variations in space and the respective maximum 
and minimum are 4.5406 10~ 20 mol./j,m 3 and 4.52854 10~ 20 mol.fim 3 which 
corresponds respectively to an increase of 3.6 tubulin-GTP heterodimers (the 
trail) and a decrease of 69 heterodimers (the depletion) in a volume of 1 /im 3 , 
if we refer to the initial concentration of 4.54 10~ 20 mol.fim 3 (equivalent to 
about 27340 tubulin heterodimers). The tubulin-GDP heterodimers produced 
at the shrinking ends are diluted in a very large amount of free tubulin het- 
erodimers. These tubulin-GDP heterodimers can be converted by nucleotide 
exchange and added to the local pool of tubu lin-GTP (rate k req at high tubulin 
concentration of about 0.02s _1 according to Melki et al (1988)). Nevertheless, 
the average increase is 3.6 heterodimers in a volume of 1 /xm 3 , that already 
contains more than 27000 heterodimers. In consequence, the differences of 
tubulin amounts surrounding the reacting ends are very small. Moreover, the 
fluctuations of concentration in the other parts of the medium (far from the 
tips) are also very low ; they are often variations of the order of fractions of 
individual molecules only, which can be also viewed as very low fluctuations 
of the probability of presence (density) of tubulin heterodimers. 
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2.2 Quantifying heterogeneity in the composition of the medium at the 
microscopic level. 

The microscopic level is the most accurate level for describing the action of 
microtubule ends on the chemical medium. Working with finite numbers of 
indivisible molecules instead of continuous - floating encoded - concentra- 
tions prevents the effects and the mistakes of interpretation as described be- 
fore. Moreover, having a realistic spatial representation of all the protagonist 
molecules present in the medium - free tubulin-GTP and tubulin-GDP het- 
crodimers and assembled heterodimers (e.g. oligomers, nuclei, microtubules) - 
allows observing spatial encumbrance effects due for example to the elongated 
shape of microtubules or to the channel formed by individual microtubules, 
which biases the direction and speed of molecular diffusion (Oddc 1998). 



Yet, VanBuren et al (20051 proposed a very accurate mechanochemical 
model of microtubule dynamics at the molecular level, but the consequences 
of these dynamics on the local molecular environment are not studied in their 
article. We are more interested here in the way tubulin dimers produced 
by the disassembling ends of microtubules diffuse around these ends and if 
they can be considered as sources of mesoscopic and anisotropic (elongated) 
heterogeneities of tubulin concentration and/or composition like the hypo- 



thetical tubulin trails that were advanced in Robert et al ( 1990| ) and Tabony 



et al ( 1999 ) to explain microtubule self-organization in the same way that well 
known trail systems like ants. 

We just mentioned that chemical trails and depletions are probably very 
weak. In particular, tubulin-GTP variations are hidden by the amounts present 
in the whole medium and, in consequence, are certainly unable to affect mi- 
crotubule dynamics. Nevertheless, local heterogeneities composed of tubulin- 
GDP heterodimers appear and could have an influence on microtubule dynam- 
ics at a more macroscopic level. We wished to quantify such heterogeneities 
and to observe their survival inversely proportional to diffusion. In or- 
der to give us the better chance to obtain tubulin trails from disassembling 
microtubules, we caricatured the best scenario and fixed a situation where 
microtubule dynamics are simplified, i.e. a situation where microtubules are 
only allowed to disassemble constantly at high rates (without interrupts in the 
disassembly such as short periods of assembly or pauses). In these simulations, 
microtubule disassembly consists in a constant release of tubulin-GDP, het- 
crodimer after heterodimer, and not in the release of tubulin coiled oligomers 
from proto-filaments nor the release of single heterodimers described by any 
stochastic based process. Assembly was not permitted. The conversion of 
tubulin-GDP into tubulin-GTP was not either permitted so as to measure 
the accumulation of tubulin-GDP over simulations. This approximation is 
however reasonable because of the very small rate of nucleotide exchange in 



tubulin (k reg = 0.02s- 1 ( |Melki et al||l988| ). 



Because they are not necessary for this study, we don't describe these 
reactions in the article. Moreover, we don't implement any microtubular me- 
chanics since at this scale level diluted microtubules are considered as rigid 
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rods due to their very high persistence length (5 mm) and because we also 
consider that microtubule-microtubule mechanical interactions are very rare 
in such diluted conditions. In this context, it is allowed to simulate only the 
diffusion of the released tubulin-GDP heterodimer from the tip of shrinking 
microtubules and conclude from the observed simulated behavior. 



2.2.1 Simulation procedure 

We use a microscopic representation of microtubule solutions where all tubulin 
heterodimers and assemblies are represented as small particles (with certain 
dimensions and volumes) suspended in an implicit solvent and move due to 
thermal agitation. Individual tubulin heterodimers are shaped as small ellip- 
soids (total dimension along the 3 axes : 8 x 6.5 x 4.6 nm). All reaction or 
diffusion events are coded in discrete time in a continuous space so as time 
steps can vary from nanoseconds to seconds and space steps can vary from 
nanoscopic levels to macroscopic levels. At each time step all tubulin het- 
erodimers and assemblies diffuse independently and we evaluate the direction 
and the length of their diffusion jumps and changes of orientations. We choose 
this approach because we do not need a good description of the molecular ki- 
netics for observing diffusion patterns. In the other case, when reactions are 
implemented, a probabilistic description of both diffusive and reactive events 
is more adapted. Although our implementation can also work in continuous 
time and semi-discrete space (not presented here) according to the Gillespie 



algorithm (Gillespie 1976 1977) and as implemented in Elf et al (2003), the 
discrete time method presented here has the advantage of being more precise 
in the description of diffusion at the microscopic level (microscopic or molec- 
ular diffusion). 

Preformed microtubules (radius r — 16.8 nm) were designed as helices with 



an angle step of 27.69 ° between two successive tubulin heterodimers (Langford 



1980). As for tubulin heterodimers they can be approximated as very elon- 
gated ellipsoids. 

Usually, in numerical models of molecular diffusion, one consider isotropic 
diffusion of spherical particles. Given the macroscopic diffusion constant of a 
particle D, one can calculate at each time step dt the jump dx of this particle 
in a random direction (between and 2tt) , using a Gaussian function centered 
on with variance V = 2Ddt. 

Nevertheless, many molecular assemblies in microtubule solutions have an 
anisotropic diffusion due to their elongated shape. For example, microtubules 
are elongated rods moving preferentially in the direction of their long axis, 
while tubulin heterodimers are quasi-spheric molecules having a quasi-isotropic 
diffusion. Following is the procedure used : 

— First, we approximate the 3D-shape of each molecule (here, tubulin het- 
erodimers or microtubule supramolecular assemblies) by an ellipsoid (Fig. 
[3l. Ellipsoids are the only quadrics that can be evaluated from sets of 3D 



points, that give information of center and orientation of an object (Li and 
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Griffiths|2004[ Han et al|2006 ), and for which the description of individual 
translational diffusion and rotational rates - expressed in terms of corre- 
lation times (i.e. the number of complete rotations during a given time) - 
along the 3 axes are described. At this spatial scale level (< 1 fim 3 ) and 
time scale level (nanoseconds to microseconds), the number of different 
molecular assemblies that form is limited and their length are also limited 
and can all be well-fitted by ellipsoids. 

Using a variation of the Stoke-Einstein formula (Einstein 1956| ), called 
the Perrin formula (Koenig 1975| |Jonstromer et al 1991), we calculate 
the translational diffusion rates (D T k) and the correlation times (D R k) 
along each axis of the considered molecule. The Stoke-Einstein formula 
Dk = kf,T/ fk is usually defined for spherical particles. Since we simulate 
the diffusion of elongated molecular assemblies like microtubules, we apply 



a shape factor Sk, called S Perrin factor (Koenig 19751, on the frictional 
coefficient f T f. — 6irr]rSk for the translation and f R k = 6rjVhSk for the 
rotation, Sk being equal to 1 for a spherical particle, and Vh is the volume 
of the equivalent hydrated spherical particle. Sk, where k designates the 
axis (x,y or z), takes into account the shape of the ellipsoid (its 3 radii) and 
its nature, i.e. oblate ((flying saucer) or prolate (cigar rod). For example, 
microtubules are like cigar rods, i.e. prolate ellipsoids, whereas tubulin 
rings correspond more to oblate ellipsoids. There are two kinds of shape 
factors : one used for balancing the translational diffusion rates (S T k), the 
other for the correlation times (S R k). Descriptions of the algorithms for 
calculating the translational diffusion rates and the correlation times along 
the three axes are given below. 

At each time step, the jump and rotation events of each molecule are 
obtained separately by two Gaussian functions along the 3 axes of the 
molecule, using the respective values of D T k and D R k for calculating the 
variances of the Gaussian curves. Then, each molecule re-orientates in the 
direction given by the calculated rotational rates, and moves according to 
its new orientation and to the calculated translational vector. Molecules 
diffuse randomly with the constraint that 2 molecules can not exist at the 
same time in the same place (their hydrated volumes can't intersect). 



2.2.2 Translational rates 

Let us consider an ellipsoid and its 3 radii r$ along its 3 axes a%, i = 1..3. The 
translational rate D T k along a given axis ak is only depending on the shape 
of the elliptic projection of the ellipsoid along this axis onto the perpendicular 
plane, as to say it depends on the values of the two radii r/, I ^ k and 
r m , m k. The approximative shape of a tubulin heterodimer is a short 
prolate ellipsoid (e.g. a cigar) such as r% — r^ > r^ — r%, with t\ > r<i > r%. 
Long and straight microtubules or linear oligomers of tubulin, for which the 
ratio p of radii is very large (e.g. p ~ 30 for a 1 /im long microtubule, about 
p ~ 8 to 10 for a tubulin ring or for a 4 heterodimers long oligomer, compared 
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to p ~ 2 for a single tubulin heterodimer) , can also be well approximated by 
long prolate ellipsoids (e.g. rods). Nevertheless, when the ratio between the 
principal radius and the secondary radius is close to 1, the projected shape is 
also not so far from that of an oblate ellipsoid (e.g. a flying-saucer), as to say 
a disk-shaped projection. This is the case for tubulin hcterodimers (ratio = 2) 
or more for tubulin rings (ratio ~ 0.36). 

Considering the axis of translation a*,, the ratio of the two perpendicular 
radii is p — r t /r m , r ; > r m , the contribution of the prolate shape (rod or 
cigar like) is S Pro = \Jp 2 — l/(p 3 ln(p + y 'p 2 — 1)) and that of the oblate 
shape (a disk) S om = ^Jp 2 — l/(p 3 arctan(\/p 2 — 1)). Ordinary shapes can 
not be always described as pure prolate or oblate ellipsoids. Then, the shape 
coefficient of any ellipsoid, S T k for a translation along Ofc, is defined here 
in first approximation as a balanced combination of both oblate and prolate 
contributions in such a way that S T k — [S OM + (p — 1) S Pro ]/p if r; > r m 
for any ellipsoid, and S T k = 1 if r; = r m in case of an exact disk-shaped 
projection. The frictional coefficient of the Stoke-Einstein formula becomes : 
f T k = 6irr]S T k^/rir m . For tubulin heterodimers or microtubules, S T i ~ 1, a\ 
being the principal axis of the molecule, because r2 ~ r%. 

2.2.3 Correlation times 

Along the principal axis a\ of an ellipsoid (the longer axis for a prolate el- 
lipsoid or the shorter axis for an oblate ellipsoid), the rotation is always of 
prolate type. The two secondary axes of tubulin heterodimers, of straight 
microtubules or of tubulin rings are equal so the shape factor of rotation 
along their principal axis cti, is S R \ — 1. To determine the type of rotation 
along the two secondary axes, we have first to determine the global shape 
of the fitting ellipsoid, as described before. Tubulin heterodimers or long 
and straight microtubules are prolate ellipsoids. Their shape factor is de- 
scribed as follows. Let's define au the principal axis of the ellipsoid and ai 
and a m the two secondary axes. If we consider the rotation along one of 
the two secondary axes, ai (resp. a m ), the ratio of its perpendicular radii is 
q = Tm/fk, ffc being always the radius collinear to the principal axis of the 
ellipsoid. The shape factor that acts of the frictional coefficient of rotation 
along the axis ai (resp. a m ) is S R i = 4(1 — <j 4 )/[3g 2 ((2 — q 2 )C — 2)] with 
C = ln((l + \f\ — q 2 )/q)/ y/l — q 2 . On the contrary, tubulin rings, for exam- 
ple, are oblate ellipsoids. The definition of their shape factor is the same than 
before with C = 2 arctan(\J\ — q 2 )/ yl — q 2 and q = r^/r m . 

2.2.4 Tubulin- GDP amount profiles 

The amount profiles are measured from the geometric center of the disas- 
sembling tips of an array of 5 microtubules (see the inset, Fig. [5] showing a 
transverse cross section of the x axis and of 5 MTs) , each of them respectively 
separated by 30 nm (one microtubule diameter) . All microtubules disassemble 
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simultaneously and regularly (no probabilistic events) at 20 \im.min~ x (about 
1.85 ms.dimer^ 1 ) which is a quite fast disassembling rate. The amount of re- 
leased tubulin-GDP heterodimers is very low (25 tubulin-GDP heterodimers 
are released by the 5 disassembling microtubules during the 11 ms of the simu- 
lation) and needs to be integrated in time for obtaining good average profiles. 
The graphic has been reconstructed by integration of the amount maps of 6 
independent simulations, during 1.8 ms (i.e. the average time separating the 
release of 2 tubulin-GDP heterodimers by a disassembling microtubule), be- 
tween the simulation times 9.2 ms and 11.045 ms (sampling time steps are 
equal to 5 10 -6 s), along the 3 axis (i.e. a total of 6642 profiles). 



2.2.5 Continuous model fitting 

To estimate the macroscopic diffusion rates from the tubulin-GDP amount 
profiles, we reproduced the numerical experiment in a continuous model of 
diffusion, in which 5 tubulin-GDP heterodimers are produced every 1.844 ms 
(5 times between and 9.22 ms) at the position x — 0. This is obtained 
by summing 5 times the equation of diffusion (Nq / (2y/n Dt)) .exp(— x 2 /(ADt)) 
from a source, each equation shifted in time by 1.844 ms. The integration 
of the signal every 5 10~ 6 s between 9.22 ms and 11.064 ms is a sum over 
this period and with the same time step of this sum of time-shifted diffusion 
equations. Two continuous models are considered : the one with only one 
diffusion parameter D corresponds to a continuous uniform diffusion model 
; the other has two diffusion parameters Dj n Bundie and DoutOf Bundle which 
correspond respectively to the diffusion of tubulin-GDP heterodimers within 
the microtubule bundle (area approximately comprised between -0.1 and 0.1 
[im] and the diffusion outside of the bundle. A large range of diffusion pa- 
rameters is explored so as to find the optimal diffusion parameters that allow 
to fit the profiles of integrated tubulin-GDP amounts. We also allowed to fit 
the profiles shifted to larger values in the limit of their standard deviation. 
The lowest sum square error between the recorded profiles and the continuous 
models indicates the optimal fit, i.e. the optimal diffusion parameters. 



3 Results 

3.1 Values of the diffusion rates. 



For individual tubulin heterodimers, we obtained the following translation 
rates and correlation times (table [2]) for water conditions as for in vitro solu- 
tions of microtubules (viscosity rjw a ter, 37 ° c = 6.915. 10~ 4 N.m.s~ 2 ) or for cy- 
toplasmic conditions as in living cells (viscosity ncytopiasm, 37 ' c = 5.6. 10 -3 N.m.s~ 2 
dSalmon et ail[l984| ). 



The same technique was applied to determine the translational diffusion 
rates and correlation times of tubulin assemblies, in particular microtubules. 
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As an example we estimated their values for a microtubule of 1 \xm long. Val- 
ues are given in table [3] From these simulations, we verified the value of 
the macroscopic diffusion of the population of tubulin heterodimers measured 
experimentally in Salmon et al (1984 1. We followed the diffusion of about 
300 tubulin-GDP heterodimers in solutions containing fixed concentrations 
of 10 mg.ml~ l of tubulin-GTP during several seconds of real time. During 
a simulation, the pathway of each molecule and its distance from its initial 
position are recorded. Toric boundary conditions are used for tubulin-GTP 
heterodimers to maintain their initial density in the sample, whereas bound- 
aries are permeable for tubulin-GDP heterodimers, to obtain their correct 
amount profiles. The estimated values obtained here were perfectly consistent 
with the values measured experimentally : D tu b, cytoplasm — 5.9 10~ 12 m 2 .s _1 
in the cytoplasm and D tubi wa ter — 48 10~ 12 m 2 .s^ 1 in water. We also car- 
ried simulations of populations of non- interacting short microtubules (1 /im 
long) in diluted solution and yield the following macroscopic diffusion con- 
stants : Dmt i fim, cytoplasm = 0.55 10~ 12 m 2 .s _1 in the cytoplasm and 

DmT 1 nm, water = 5.6 10 -12 771' 



3.2 Formation of tubulin-GDP heterogeneities in the medium. 



To test the formation of tubulin trails, we began by positioning a single mi- 
crotubule in the middle of a simulated sample of cubic shape of 1.4 /im side, 
oriented along the X axis, in a medium containing 10 fiM tubulin-GTP (5400 
heterodimers. /im" 3 J] The microtubule was not allowed to diffuse and stayed 
located on the x axis. In a first simulation, the microtubule was disassembling 
100 x faster than a real microtubule (normally shrinking at a maximum rate 
of about 20 /jm.min -1 ), as to say at a rate equivalent to 2000 ^m.min -1 . We 
realized this numerical experiment in unrealistic conditions so as to obtain 
approximately what was expected in Tabony et al ( |1999 ); Glade et al (20021 



Tabony (2006); Robert et al (1990). In this case indeed, an eye-observable 
concentrated area formed around the shrinking end (Fig. [2]). At realis- 
tic rates of disassembly (20 ^m.min -1 ) however, the variation of tubulin-GDP 
amount around the tip was very hard to detect, in particular in water. Indeed, 
molecular diffusion is a very fast homogenizing process compared to the micro- 
tubule disassembly process that creates heterogeneities. During the release of 
one tubulin heterodimer in the cytoplasm by a microtubule, the last released 
tubulin-GDP heterodimer has sufficient time (about 1.85 ms) to explore an 
average sphere of 150 nm of radius (0.5 /jm in water). In consequence, during 
microtubule disassembly - or assembly (Oddc 1997) - the solution is rapidly 
homogeneous at the micrometer scale (< 1 — 2 /im). Nevertheless, although 



2 Simulations have shown that tubulin concentrations of 100 /iM or the presence of nu- 
merous microtubules distributed in the sample doesn't change the macroscopic diffusion 
rate of free tubulin or small molecules (results not shown) . |Liu et ai] ( |2006} also mention 
that the packing geometry of microtubular bundles doesn't avoid tubulin diffusion within 
the bundles. However, our simulations show that increasing concentrations of microtubules 
modify locally the diffusion rates of tubulin heterodimers. 
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it's not intense, the very weak gradient of tubulin-GDP heterodimers can be 
measured in the simulations (Fig. [5J. The gradient shape obtained after in- 
tegration of numerous simulations is close to a sum of Gaussian distributions 
(solutions of the diffusion equation that started at different times, i.e. when 
tubulin-GDP heterodimers are released by microtubules) and roughly speak- 
ing, the tubulin-GDP heterodimers are nearest to the tip than very far away. 

This can be better observed by allowing several microtubules to disas- 
semble in the same place. We simulated 5 aligned and parallel microtubules 
disassembling at 20 /im.min^ 1 during 11 ms. A simulation time of 11 ms is 
a compromise so as to have a sufficient amount of tubulin-GDP released by 
the 5 disassembling microtubules and to have a limited diffusion area (less 
than 2 m) to avoid simulating huge volumes. This time, the formation of a 
gradient of tubulin-GDP was eye-observable and measurable. After 11 ms of 
reaction in cytoplasmic conditions, tubulin-GDP heterodimers are observed 
at a distance of about 0.7 /jm from the tips of the microtubules from where 
they were released and the gradient has an average maximum value of 4.5±1.8 
heterodimers at the disassembling microtubule tips (Fig. [5] A). In water, the 
intensity of this gradient is weaker. After 11 ms, the gradient of tubulin-GDP 
extends more rapidly until a distance (observable tubulin-GDP heterodimers) 
of about 1.5 fiui from the tips of the microtubules. It has an average max- 
imum value of 1.31±1.08 heterodimers at the tips of the microtubules (Fig. 
[5j B). As seen in table |1J the values of diffusion estimated from the aver- 
age recorded pathways of tubulin-GDP, from the uniform diffusion continuous 
model and microtubule dependent diffusion continuous model are compara- 
ble. However, the values given in table [4] and the model fits shown in Fig. 
[5] clearly confirm that only microtubule concentration dependent models can 
describe correctly the tubulin-GDP amount profiles shown in Fig. [5j This 
result disagrees with Liu et al ( 2006[ ) : the packing geometry of microtubular 
bundles seems tp avoid tubulin diffusion within the bundles. Microtubules, 
by their assembling or disassembling activity can produce local chemical het- 
erogeneities. Unfortunately, even if dozens of disassembling microtubules are 
located in the same place, the increase of total tubulin (tubulin-GTP and 
tubulin-GDP) concentration is undetectable: the variation of tubulin-GDP 
amount is indeed much weaker than the natural fluctuations of total tubulin 
amount (there are about 1000 tubulin-GTP heterodimers per released tubulin- 
GDP hctcrodimer). These heterogeneities are not concentrated areas but more 
areas in which the composition is changed. 

For individual microtubules, the heterogeneities are very weak and ex- 
tended in space. In consequence, particularly in in vitro solutions (viscosity 
of water), it is very improbable that an individual microtubule can influence 
another one in its neighborhood. The effect can become significant when 
produced by a group of synchronously reacting microtubules. To produce an 
intense composition or concentration heterogeneity, there are two possible sce- 
narios: (1) in a solution of randomly distributed microtubules, several dense 
nodes of microtubule ends naturally exist that can form initial local composi- 
tion (or concentration) heterogeneities; (2) microtubule arrays (bundles) form 
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locally by another mechanism (for example by static interactions (Liu et al 
Baulin et al||2007 >). Then, in both cases, the formation of heterogeneity 



2006 



nodes can potentially provoke the nucleation, the recruitment or the inhibition 



Glade et al 


(2002 


); 


Robert et al 


(1990) 



over, because of the rapid diffusion of tubulin heterodimers in comparison to 
the reactivity of microtubules, tubulin 'trails' are not directional as proposed 
in the reaction-diffusion ant-based model presented in Tabony et al (19991; 



Glade et al 


(2002 


); 


Tabony 


(2006) 



As their shapes seem isotropic, the trails 
produced by microtubules might be inefficient to serve as guiding chemical 
rails for other growing microtubules. 



Other effects could reinforce a little the intensity of tubulin-GDP concen- 
trated areas. First, we observed a effect of microtubule density on local macro- 
scopic diffusion : diffusion of tubulin heterodimers is about 3 times slower 
within bundles. Microtubule bundles should then contain concentrated stocks 
of tubulin-GDP, in particular in bigger bundles. Moreover, microtubules can 
release individual tubulin-GDP dimers but also oligomers of several assembled 
tubulin heterodimers. As calculated by our algorithm, their diffusion rate is 
about 2 times slower than that of free tubulin (for example, a straight oligomer 
of 10 subunits have the following diffusion rates D oUgo w subunits, cytoplasm = 
2.8 10" 12 rn 2 .s _1 and D oHgo w subunits, water = 23 10~ 12 rn 2 .s _1 compared to 
Dtub, cytoplasm = 5.9 10~ 12 m 2 .s^ 1 and D tu b, water = 48 10~ 12 m 2 .s _1 ). This 
could help a little to maintain a little the free tubulin-GDP more concentrated 
in the neighborhood of the tip. 



3.3 Effect of microtubule diffusion. 

In the simulations described in the previous chapter, microtubules were as- 
sumed to be motionless (cf. Fig. [5]). In real solutions, they can diffuse 
when sufficiently small, not blocked by other interacting microtubules and 
not bounded to the walls of the sample. Their diffusion is rapidly limited 
by their size, but the diffusion of microtubules of about 0.1 to 1 fim long or 
more is not negligible. While tubulin particles have an isotropic diffusion, in 



agreement with Han et al (2006), microtubules have a 'short time anisotropic' 
diffusion and a 'long-time isotropic diffusion' due to their anisotropy of shape 
(see Fig. [6] showing records of individual trajectories of tubulin dimers or 
of microtubules). Their diffusive motion will then spread more the patterns 
of tubulin-GDP diffusion^] obtained in Fig. [5] by changing the places where 
tubulin-GDP is released (Fig. [7]). We did not estimate the macroscopic diffu- 
sion rates in this condition because we would need a much more complicated 
continuous model for, in which the microtubule dependent diffusion area would 
change over time. In these simulations, microtubules were not grouped in 



3 The diffusion constant that describes the spreading of this pattern (the trail) results 
from coupled tubulin-GDP and microtubule diffusion but also depends on microtubule dis- 
assembling rate. 
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the form of bundles. They could interact and collide but had no cohesion. In 
real solutions, microtubules can form bundles densely populated where micro- 
tubules stay grouped, behaving more like in Fig. [5j As seen before, this may 
reduce this spreading effect due to microtubular diffusion. 



4 Discussion 



In this article, we asked the question of whether two (or more) microtubules 
close together in a solution of tubulin can communicate by way of chemical 
interactions as observed in collective systems like ant colonies and whether this 
can explain microtubule spatial self-organization from microscopic to macro- 
scopic levels. 

Foraging ants normally walk randomly but they can let behind them con- 
centrated trails of chemicals called pheromones. This communication way is a 
manner by which is recorded the pathway of ants between a source of food and 
their nest once food has been found. When neighboring ants cross these trails, 
the probability they use the signalized pathway is greater than continuing to 
walk randomly. This communication way is very efficient because the chemi- 
cals released by these insects are very concentrated locally, even after a long 
period of diffusion (and degradation) , because the dispersion area by diffusion 
around the original trail is of the same order as that of the size of the agents 
(the ants), because a trail of pheromones is an isolated signal, not so much 
diluted in an environment full of equivalent pheromones, and because, due to 
these reasons and because social insects are sensible to weak concentrations of 
pheromones, neighboring ants use efficiently these trails, release pheromones 
within the same trail, and the signal reinforces. 

On the contrary, the present numerical simulations have shown that the 
tubulin trails are of very low intensity, that the magnitude of diffusion of these 
'trails' is largely greater than the size of microtubule ends, that the signal of 
the trails is completely diluted in an environment full of tubulin heterodimers, 
and that in diluted solutions microtubules continue to diffuse randomly, don't 
stay growing within the trails and do not work at reinforcing them. This is 
what our numerical simulations predict until the tubulin trails will possibly 
be observed experimentally or not. Moreover, even if dozens of microtubules 
produce a common trail, this would be certainly too weak to change locally - 
at a microscopic scale - the reactivity of neighboring individual microtubules. 
It results that the explanation based on tubulin concentrated trails that con- 
stitutes the basis of microtubule self-organization at the individual level in 
stationary self-organized microtubule solutions ( |Robert et al (1990); Tabony 



et al (1999); Glade et al (2002); Tabony (2006)) is probably unfit. 



Self-organization as observed in Tabony and Job (1992b 1990); Tabony 



(J1994J); |Papaseit et al| (|1999j); |Tabony et al| ( |1999[ ) is, of course, a kind of 
dissipative process : one can not deny that the reactivity of microtubules is 
implied during the development of the self-organizing pattern and its main- 
taining. However the description of that system is incomplete because it only 
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refers to a reaction-diffusion system expressed in terms of Turing-type pro- 



(Robert et al 1990 


Tabony et al 1999 


Glade et al 2002 


Tabony 2006) 



microtubule solutions, there are reactive processes and transport of matter, at 
least by molecular diffusion processes, but other processes such as biomechan- 
ical or electrostatic interactions should also have very strong contributions. 
The question now is to know how strong the coupling of reaction and diffusion 
processes is compared to the influence of static ones. On this point, our sim- 
ulations can not yet answer. It is indeed plausible that static processes play 
the most important role in the microtubule self-organizing process and that 
the existing tubulin-GDP trails (estimated in this study) reinforce it or even 
biased it to cause cause complex morphologies. 

In their paper, Liu et al ( 2006 ) analyzed how the microtubular self-organized 
stripes formed. They suggest that microtubules are packed into bundles and 
buckle due to their growth. In this scheme, the self-organized microtubule 
stripes form by this mechanism from previously aligned microtubule bundles 
that are initiated either by static magnetic fields or convective flow (shearing). 
This mechanism is similar to those proposed by Portet et al (2003); Baulin 



(2003); Ziebert and Zimmcrmann (20041; Baulin et al (20071, that couple mi- 
crotubule growth to microscopic self-ordering due to nematic ordering, the 
macroscopic self-ordering being biased by gravity. 

A kinetic experiment of neutron scattering on a microtubular self-organizing 
structure examined through a horizontal slit (dimensions : 4 x 0.5 mm), de- 
scribed in Tabony (1994) (see Fig. 4 of this article) and obtained in mi- 



crotubule stationary self-organized morphologies, suggested that microtubules 
disassemble and reassemble during the organizing process. Because of this, mi- 
crotubules are though reorient and form progressively the patterns of concen- 
tration and orientation. Microtubules were oriented in a preferential direction 
within the first 2 hours, then the time dependence of microtubular scattering 
intensity showed a decrease at about 2 hours, indicating a loss of any preferen- 
tial orientation, and then increased again indicating a reorientation along the 
opposite direction. However, despite of the interpretation given by the author, 
it does not indicates inevitably a 'partial disassembly-reassembly-disassembly' 
process, nor it gives information on what really occurs at the microscopic level. 
This interpretation would be right granted that microtubule packs can't move 
in the solution. It only indicates that microtubule orientation changes and 
that, at one moment, their orientation seems isotropic in the field of obser- 
vation. This result could be interpreted differently, in agreement with the 
bundling and nest buckling process mentioned before : packs or single bun- 
dles of aligned microtubules buckle due to microtubule growth. This causes 
changes of the orientation of the bundled microtubules and also causes their 
mechanically-driven travel (in opposition to the treadmilling travel of micro- 
tubules due to their reactivity) across the solution and the observation window. 



At the beginning of the article, we mentioned that different mechanisms could 
be involved in the different varieties of microtubular self-organizing behav- 
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iors. On the one hand we have chemical solutions in which microtubules are 
very dynamic (as observed in time series of the concentration of assembled 
tubulin measured by spectrophotometry at 350 nm) and that clearly imply 
reaction coupled to diffusion processes ensuring synchronizing phenomena in 



Carlier et al 1987 Mandelkow 



microtubule populations (Pirollet et al 1987 
et al||1989 1 - although it does not imply that static processes of ordering are 



absent. On the other hand, spatially self-organizing solutions of microtubules 



(Tabony and Job||1990 Hitt et al|1990 


Tabony and Job||1992b 


Tabony|1994 


Papaseit et al||1999 


Tabony et al||1999 


Liu et al||2006 ) are less dynamic (only 



one overshoot of microtubule concentration is observed) and probably involve 
other contributing mechanisms such as biomechanical and/or nematic order- 



ing (] 


Hitt et al|1990 


Portet et al|2003 Baulin||2003 ; Zicbert and Zimmermann 


2004 


Liu et al 


200( 


3). However, the reactants are very similar from one of 



these 

concentration of magnesium ions (20 mM in Mandelkow solutions, about 10 
mM in those of Carlier, Pirollet and Hitt ones, and only 1 mM in Tabony and 
Hitt solutions), the use of deuterium oxide instead of water only in Tabony 
solutions, and the presence of stabilizing MAPs in Hitt solutions. Magnesium 
ions are known to increase the reactivity of microtubules : they promote mi- 
crotubule assembly and disassembly ( Flyvbjerg et al||1996 ). On the contrary, 
deuterium oxide has been shown to stabilize microtubules^ ( Chakrabarti et al 
1999 ) : it suppresses the dynamic instability and the treadmilling behavior of 
microtubules ( Panda et al||2000 | but stimulates the nucleation of new micro- 
tubules from free tubulin heterodimers ( Itoh and Satofl98"4 ). In microtubule 
solutions where stationary self-organized morphologies appear (Tabony and 



Job 1992b) everything tends to produce and maintain numerous, very stable 
and probably very long microtubules whereas it is exactly the contrary in the 
case of very reactive solutions that favor the formation of temporal oscillations 
of of traveling waves of microtubule concentration. 

This re-examination of experiments on microtubules over the last 20 years 
gives arguments, which add to the numerical results presented here, against a 
'molecular ants'-based scheme and against a reaction-diffusion scheme for the 
stationary well organized spatial morphologies desc ribed in|Tabony and Job 



(1990); Hitt et al (1990); Tabony and Job (1992b I; Tabony (1994); Papaseit 



et al (1999); Tabony et al (1999); Liu et al 



(20061. On the contrary, in this 



case, the thesis of a self-organizing mechanism based on the biomechanics of 
individual growing microtubules and bundles appears more reasonable. 

Nevertheless, in microtubule solutions showing periodic temporal varia- 
tions or spatio-temporal variations (propagating waves) of the concentration 
of microtubules, the contribution of tubulin-GDP 'clouds' produced by indi- 
vidual microtubules or by bundles of numerous microtubules can't be ignored. 
In these solutions, microtubules are synchronized within the bulk solution gen- 
erating temporal oscillations of the microtubular concentration of the whole 



4 D2O also prevents tubulin proteins against denaturation, this allowing to realize long 
term experiments 
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solution or within distances of the order of 1 mm allowing to form propagat- 
ing waves. This may be mediated by millimeter scale variations in the tubulin 
composition of the medium. The range of the clouds of tubulin-GDP pro- 
duced by dense nodes of microtubules should control the distance at which 
the reactivity of 'neighboring' microtubules can be modified and the micro- 
tubules to be synchronized. The present study only reports how tubulin-GDP 
heterodimers diffuse from simplified microtubules. Further advances including 



reactions based on rescue-catastrophe dynamics (Surrey et al 2001 Nedelec 



2002 Nedelec et al 



2003) or more accurate kinetics such as those described 



VanBuren et al (20051 will allow understanding, at the microscopic level, 



what controls microtubule synchronization over mesoscopic to macroscopic 
distances. 



In this article, further than the study realized on microtubule self-organization, 
we addressed the question of the robustness of a microscopic process (i.e. here, 
the theoretical communication between microtubules by way of tubulin trails) 
that, confronted to ambient noise or accidents, could drive the elements of 
the system to the appearance of macroscopic order. In biology, physics or 
biocomputing, more and more studies deal now with this question where the 
dynamical self-organization of a system depends on collective dynamics based 
on microscopic processes of communication, often in a noisy context (i.e. ther- 
mal agitation, external perturbations ...) (Lcsnc 2008). If Nature works since 
billions years to produce and select processes that are robust in such perturbed 
contexts for realizing a coherent - macroscopic - function that plays a role in 
living organisms, this is not the case for our engineered products. Now, theo- 
reticians and engineers expect from nano or microscopic natural - biological, 
physical - systems or from systems inspired from nature (e.g. synthetic fiber- 
shape assemblies Rothemund et al| (2004a b); Glade (20081) to realize tasks 
or computations. This knowledge presents a real interest in collision based 
computing ( Costello and Adamatzky 2005 ; Igaras hi et al|2006 ) or in dynamic 



self-organizing molecular processors (Pfaffmann and Conrad 



2000 Tcuscher 



2007 Glade et al 2009) in which the wires are not as clearly defined as in 



electronic processors or brains (Demarse et al 2001) but constituted by the 
agents themselves and by their pathways ( Glade et al|2009 ). The spatial scale 
level is certainly the most critical criterion that controls the efficiency of these 
self-ordered microscopic processes to realize the expected macroscopic actions 



( Conrad 1995 ) or to generate the expected processing architectures ( Pfaffmann 
and Conrad||2000| |Teuscher|[2007| |Glade et al||2009[ ). Studies of how in forma- 



tion or matter is exchanged between agents (as presented here or in Lizier 



et al ( 2008 ) ; Glade et al ( 2009 1 of similar systems are a necessary preliminary 



to such unconventional computing or robotic works, but also concern strongly 
all the agent-based models of molecular and supramolecular organization in 
cellular biology, a field that tends to replace (or at least complement) increas- 
ingly the differential equation based models. Our work, although not in favor 
of a reaction-diffusion based computation with microtubules, could be used in 
such a way, for determining the efficiency of information or matter exchange 
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in similar systems such as actin comets or other potential trail systems, i.e. 
chemotactic cells ( Bagorda et al|[2006 ). 
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Table 1 Experimental conditions that allow microtubules to self-organize spa- 
tiotemporally, temporally, or spatially. 
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Fig. 1 Hypothetic self-organization scheme at the level of individual micro- 
tubules as proposed in Gla de et al| ( [2002[ l: |Tabony| ( |2006[ ). In this scenario, micro- 
tubule disassembly generates concentrated trails of tubulin-GDP (in black) rapidly regen- 
erated into tubulin-GTP (in gray). These concentrated trails are favorable media for the 
preferential growth (or nucleation) of neighboring (or new) microtubules. The consequence 
is the formation of microtubule arrays. 
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Fig. 2 Tubulin trail in a reaction-diffusion model. This tubulin trail is simulated 
using a 2D spatiotemporal differential equation system modeling a reaction-diffusion micro- 
tubule system (dimensions: 200 X 100 fim) similar to that proposed by |Glade et al| ( |2002[ ) . 
The figure shows the tubulin-GTP concentration map (represented in false 3D and grey 
scales), in superposition of a single microtubule (above). This microtubule measures 25 fxia. 
In this simulation, an isolated microtubule is in treadmilling motion, growing at the right 
end (at 21 fim.min -1 ) and shrinking at the left one (at 18 fim.min~ ), in a medium con- 
taining initially 5 mg.ml^ 1 - or 4.54 10~ 28 raol.(im 3 - of tubulin-GTP. At the growing end, 
an 'intense' depletion of tubulin forms (the black hole). At the shrinking end, tubulin-GDP 
is released and rapidly converted into tubulin-GTP, producing a trail (white trail at the 
left of the depletion). The hole of the depletion is more intense than the peak of the trail 
because tubulin-GDP heterodimers have to be regenerated first into tubulin-GTP at a rate 
of 0.2 s~ . During that time they diffuse all around lowering the intensity of the trail. 




Fig. 3 Approximation of molecular assemblies by ellipsoids. (Top) an approxi- 
mation of rod-like microtubules (left) and quasi-spherical tubulin heterodimers (right) by 
ellipsoids is shown. (Down) This allows correcting the translational rates of diffusion and 
their correlation times of the particles in solution by a shape factor calculated along each 
axis. 
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Table 2 Diffusion rates and correlation times of tubulin heterodimers. These dif- 
fusion parameters calculated for tubulin heterodimer elliptic approximate take into account 
the viscosity of water or of the cytoplasm at 37 ° G. 
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Table 3 Diffusion rates and correlation times of a microtubule of 1 fim long. 

These diffusion parameters calculated for the elliptic approximate of a 1 /im long microtubule 
take into account the viscosity of water or of the cytoplasm at 37 ° C. The microtubule is a 
cylinder so its axes y and z are equivalent. 
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Fig. 4 Concentrated area of tubulin-GDP forming at the shrinking end of a 
fast disassembling microtubule. Tubulin-GTP heterodimers are displayed in gray while 
Tubulin-GDP is in black. In this simulation, the microtubule disassembles 100 x faster (at 
2000 /im.min^ 1 ) than a normal microtubule. 
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A B 

Fig. 5 Profiles of tubulin-GDP amounts around the tips of 5 disassembling 
microtubules. (A) Profiles produced by 5 motionless microtubules (the inset shows a 
transverse cross section of the x axis and of 5 MTs), in cytoplasmic conditions. The macro- 
scopic diffusion rate of individual tubulin hctcrodimers measured fr om their trajectorie s 
corresponds to that measured in the cytoplasm (5.9 10 — 12 m 2 .s — 1 ) jSalmon et al|[l984| . 
(B) The same in water at 37 ° C. The average value of the macroscopic diffusion constant 
measured from their trajectories (4.9 10 -11 m 2 .s _1 ) is about 8 times larger than in the 
cytoplasm, as measured in water by |Salmon et al] ( |1984| | . The resulting amount profile of 
tubulin is weaker but detectable with 5 disassembling microtubules. In both graphics, the 
fits produced by continuous models based on a uniform diffusion or on a microtubule de- 
pendent diffusion are drawn. The best fit of the profile shapes can be obtained after shifting 
the profiles to bigger values within the standard deviation of the profiles. 
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Table 4 Macroscopic diffusion estimated from the profiles shown in figure|5| Dif- 
fusion rates are estimated from the average recorded pathways of tubulin-GDP heterodimers, 
or from continuous models based on a uniform diffusion or a microtubule dependent diffu- 
sion. In the latter, diffusion of tubulin-GDP heterodimers is different inside and outside the 
microtubule bundle. Sum square errors (SSE) between the model fit and the profiles are 
indicated. 
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(SSE : 0.015) 
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Fig. 6 Time-record trajectories of microtubules and tubulin heterodimers. Pro- 
jections, from the 3D space to the YZ (A) and XY (B) planes, of 3D trajectories of 3 
individual tubulin heterodimers (orange, blue and violet) and 2 microtubules of 0.1 fim 
(green) and 1 am (red). Microtubules are initially oriented along the Z axis (perpendicu- 
larly to the plane XY). At the left bottom of each trajectory record, a scheme indicates the 
initial orientation of microtubules. The trajectories are relative to their starting point and 
were recorded during 11 ms (B) and 3 s (B). Sampling time steps are equal to 5 10 — 6 s. 
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Fig. 7 Profiles of tubulin-GDP amounts released by an array of 5 disassembling 
and diffusing microtubules. Microtubules (of initial length equal to 0.5 fim long) are 
allowed to diffuse in water (A) or in the cytoplasm (B) during the simulation (11 ms). Both 
are measured as described in fig[5]from the initial geometric center of the disassembling tips 
of microtubules (denned at the beginning of the simulation). 



